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ABSTRACT 

We study dynamical effects of cosmics rays (CRs) on the thermal instability in the 
linear regime. CRs and the thermal plasma are treated as two different interacting 
fluids, in which CRs can diffuse along the magnetic field lines. We show that growth 
rate of the magnetothermal condensation mode is reduced because of the existence of 
CRs and this stabilizing effect depends on the diffusion coefficient and the ratio of the 
CRs pressure to gas pressure. Thus, a slower rate of structure formation via thermal 
instability is predicted when CRs are considered. 
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1 INTRODUCTION 

Birtli of the stars and the complex observable patterns 
within the interstellar medium (ISM) are mainly understood 
by the instabilities that may start with small amplitude per- 
turbations le ading to highly nonlinear processes, i.e. tur- 
bulen ce (e.g.. Sc alo fc Elmegreenll2004l : iMac Low fc KlessenI 
|2004| ). In fact, among the most important physical ingredi- 
ents in any theory of star formation one can mention self- 
gravity, magnetic field and net cooling function. However, 
' there is another significant physical factor that its dynam- 
ical role has not been considered in most scenarios of star 
formation until rather recently: Cosmic Rays (CRs). 

Equipartition exists between energy densi ties of thermal 
, gas, magnetic field, turbulence and CRs (e.g.. lGaisserlll990l : 
. lFerrierg|200ll ). The energy equipartition encourages one to 
' expect significant dynamical effects of CRs on the structure 
formation i n ISM as has been realized by Parker in a simpli- 
■ fied model (|Parkeil 19661 ). It is believed that giant cloud com- 
plexes in the sp iral arms of galaxies are forming via Parker 
instability (e.g.. lMouschoviaall974l : lElmegreen fc ElmegreenI 
Ll986). More recent studies show that the growth rate of the 
Parker instability is very sensi tive to the CR pressure and 



the diffusion coefficient 



Hanasz 1997: Ha nasz fc Leschllloooi n Kuwabara et al.l 



Kuznetsov fc Ptuskin 



198a 



2004 



kuwabara fc Kal2"006l ). CRs have also a vital role in propa- 



gation of MHD waves. ILo fc Kol (|2007l ) studied the the sta- 
bility of a cosmic-ray plasma system. This stability analysis 
has been generalized to a four-fluid cosmic-ray-MHD system 
which comprises magnetized thermal plasma, cos mic rays, 
forwa rd and backward propagating Alfvn waves (|Ko fc Lol 
l2009l ). 
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Interactions of CRs may al so operate as a global heat- 
ng m ec hanism in ISM (e.g.. [ Field. Goldsmith fc Habing 



19691: iGoldsmith. Habing fc Fieldl Il969l ). Recently, 
Yusef-Zadeh. Wardle fc Rovl (|2007l ) proposed that due 



to heating of the gas clouds in the central regions of our 
Galaxy by an enhanced flux of CR electrons, the rate of 
star formation in these star forming sites decreases and a 
similar mechanism may operate i n the nuclei of the other 
galaxies. In another related study, iFatuzzo. Adams fc Melial 
(2006) studied the interplay between molecular clouds and 
the irradiation by CRs produced by supernova remnants. 
This increased CRs flux has important consequences for 
star formation. In particular, a higher ionization level due 
to CRs lead to a longer ambipolar diffusi on time scale an d 
therefore slower star formation rate (Fa tuzzo et al.l [20061 ). 
However, none of these studies considered the possible 
dynamical effects of CRs on the star forming regions. 

A linear analysis of non-magnetized thermal gas 
and CR has revealed that CRs suppress the growth 
of small amplitude pertur b ations jKuznetsov fc PtuskinI 



19831: iBegelman fc Zweibel Il994l : IWagner et al.l |2005| ) 



Kuznetsov fc PtuskinI l|l983D (hereafter KP83) studied sta- 
bility of ISM consisting of regular and turbulent magnetic 
flelds, thermal gas and CRs. Although stabilizing effect due 
to the existence of CRs has been shown, heating-cooling of 
the system has not been considered by KP83 and a poly- 
tropic equation of state is used instead of the energy equa- 
tion for the gas component. Moreover, their analysis resem- 
bles to the classical Parker instability, in which the initial 
unperturbed state is determined by the balance of the forces 
rather than thermal equilibrium stat es as we will consider 
in our study. Therefore, analysis of iKuznetsov fc PtuskinI 
1 19831 ) can not address stability of a system with CRs when 
the net cooling function is playing a signiflcant role like 
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structure forma tion in HI reRions due t o the thermal con- 
densation (e. 
IWagner et al 



.Fukue fc Kamavall2007f ). In another study, 
derived stabihty condition of a non- 
magnetized thermal system for a few simplified cases. How- 
ev er, growth rate of the perturbations has not been studied 
bv lWagner et al] l|2005l ). In this study, we follow a similar 
approach but in a magnetized case considering the net cool- 
ing function of the system and the dynamical effects of CRs. 
In this regard, our study is different from KP83 who did not 
take into account thermal behavior of the system and as we 
will show growth rate of the perturbations and the criteria 
of stability are completely different in comparison to KP83 . 
Also, our study is complementary to IWagner et alj (|2005l ) 
who did not consider magnetic field and the growth rate 
of the perturbations has not been calculated in detail. It is 
assumed that CRs particles can diffuse along the magnetic 
field lines (see also. ISchlickeiser fc Lerch3ll985l ). In the next 
section basic equations of a magnetothermal system includ- 
ing CRs are presented. A dispersion relation is obtained for 
the linear perturbations in section [31 Analysis of the the 
growth rate of the condensation modes is done in section 
|4l We conclude by a summary of the results and possible 
implications in the final section. 



2 GENERAL FORMULATION 

There are three different approaches to study the dynam- 
ics of CRs. In the particle-particle approach, the plasma 
and CRs are considered as particles that may interact with 
each other via complicated processes. In a simpler approach, 
known as fluid-particle, the plasma is treated as a fluid, 
though CRs are still described as particles. The simplest 
approach is the fluid-fluid approach in which CRs and the 
thermal gas are described by different interacting fluids. The 
hydrodynamic approach can not provide the spectrum of 
CRs, however, it is a good approxima tion to use to analyze 
dynamics o f a plasma with CRs (e.g., iDrurv fc Voelk|[l98ll : 
iDrurvl 19831 ). We adopt the hydrodynamic approach to study 
effects of CRs on the thermal instability; 

We also follow the same steps as in lFieldl l| 19651 ) , except 
that CRs are considered as a fluid, and diffusion is consid- 
ered only along magnetic field lines. For simplicity we neglect 
ionization and heating by CRs, since their effects in the ab- 
sence of the dynamical role of CRs are well understood (e.g., 
[Goldsmith et al.lll969l ). 

The basic equations are 

dp „ 

It = 

P^ = ^(VxB)xB-V(p + p„), 

dB 

'dt 

V.B = 



(1) 

(2) 

(3) 
(4) 

1 ^^^^^^ 
J — 1 at 7 — 1 p di 

and 



= V X (v X B) 



where d/dt — d/dt + v.'V is the Lagrangian time derivative, 

r = -K||b(b.Vpcr), (7) 

is the diffusive fiux of cosmic-ray energy and K|| is diffusion 
coefficient along magnetic field lines. All the variables have 
their usual meaning. Also, b is a unit vector along the mag- 
netic field lines, i.e. b — B/_B. The adiabatic indices of the 
thermal gas and cosmic rays are denoted by 7 and 7cr, re- 
spectively. Also, Q represents the energy losses minus energy 
gains per unit mass. The coefficient of thermal conductivity 
K has the values K\\ and K± in directions parallel to and 
perpendicular to the magnetic field B. Finally, we can write 
equation of state as 



R 
—I 



-pT, 



(8) 



1 dpcr 7cr Per dp 



7cr — 1 dt 7cr — 1 P rff 



v.r = 0, 



(6) 



where R is the gas constant and fi represents the molecular 
weight. 

The energy equation of CRs in its complete form in- 
cludes also an extra term correspo nding to an effective 
CRs e nergy loss (see equation (13) in ISchlickeiser fc Lerchd 
(1985')). But this term is neglected in our analysis simply 
based on a time scale argument. We neglect the ionization 
energy losses of the cosmic rays because these are surely on 
a much longer time scale than we are interested. In fact, 
CR pressure is mainly from mildly relativistic protons, and 
that these have very low ionization energy loss rates and 
also low nuclear collision rates. Therefore, typically the ion- 
ization and nuclear loss time scales are around 10* or more 
years for normal interstellar values. But typical cooling time 
scale is around 10® years. Also, there is an effect due to 
second-order Fermi ac celeration of cosmic r a ys ( s ee equa- 
tions (15) and (13) in ISchlickeiser fc Lerchel |l985l )). which 
in fact would turn the cosmic ray loss time into a cosmic 
ray gain time. This effect is not so easy to discard, as this 
acceleration is caused by the same MHD waves that con- 
trol the spatial di ffusion of cosmic rays. In fact it is known 
l|Schlickeiseilll989l ) that the product of spatial diffusion time 
To and the acceleration time scale rp equals a constant that 
depends on the spectral index of the turbulence spectrum 
and (L/va)'^, where L denotes the size of the system and 
is the Alfven speed (see equation (86) in lSchlickeiseij (|l989l )). 
The time scales ar e related approximately totf = (L/va)^ 
(jSchlickeisej Il989l ) , where the spatial diffusion time scale 
is td ~ L^/ii. For a typical cooling ISM we can assume 
L = 10^" cm. Bo = lO"'' G and the ion density is n, = 10"^ 
cm~^. So, we obtain td = 10^^ s and tf = 7.2 x lO^'' s 
and so, rp S> td which means we are indeed in a parameter 
range to neglect the term corresponding to an effective CRs 
energy loss. Thus, it seems our approach is fine. We actually 
have no source term either, but we are interested in the effect 
of a pre-existing and stable CR population on the magne- 
tothermal condensation modes. To introduce additional time 
scales related to sources and sinks is a nonsense unless they 
are comparable to the time scales we are studying. 

Note that in the hydrodynamic approach of CR, the 
energy spectrum of the particles of CR is not considered in 
detail. Moreover, type of the CR particles is not specified 
in this analysis and they can consist of electrons or protons 
in any combination. In comp arison to th e classical magne- 
tothermal instability analysis (|Fieldlll965h . there is an extra 
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pressure term due to the CRs in the equation of motion ([2]) . 
So, one may expect an stabilization effect of CRs. But equa- 
tion ^ describes advection of CRs by flowing gas and the 
diffusion of CRs along the magnetic field lines through the 
thermal gas. Little is known about the diffusion coefficient 
K|l and its possible depend ence on the physical variabl es of 
the system. As explained in lSchlickeiser fc Lerchd (|l985l ) the 
parallel spatial cosmic ray diffusion coefficient k\\ in a cos- 
mic ray fiuid theory is an effective diffusion coefficient aver- 
aged over all cosmic ray momenta. Depending on the actual 
momentum distribution function of cosmic rays and the mo- 
mentum dependent cosmic ray diffusion coefficient, its value 
can be ve ry different from the standard value 10'^^ cm^ s~^ 
(see also, teerezinskn et al.iri990l : |Ptuskinll200ll '). 



3 LINEAR PERTURBATIONS 

We completely neglecting gradients in the background 
medium and the velocity is zero in equilibrium state po, 
To, Po, PcrQ. Also, we assume Q,{po,To) — 0. By perturbing 
of the form 5X{r,t) = (5Xexp(aji -I- ik ■ r), the equations 
([T])-® become 



ujSp + jpok.<5v = 0, 



(9) 



a;po5v + ik5p + ik5pcr + i(Bo.5B)^ -i(k.Bo)^ = 0, (10) 

47r 47r 



ujS'B + iBo(k.5v) - i(k.Bo)5v = 0, 



-Sp ^ 



7-1 (7 - l)po 



Sp + poVlpSp + poVItST 



(11) 



(12) 



LU I^7crPcr0 c- , (k.Bo)^ , , 

"''Per ^ 77 lT—"P+'^\\ OPcr = 0, (13) 



(7cr - 1)P0 



5p 5p ST 







(14) 



po po To 

Note that the derivative flp = {dQ/dp)T and fir = 
{dil/dT)p are evaluated for the equilibrium state. 

Then, we introduce the coordinate system e^, Cy, and 
fiz specified by 



Bo 



Bo X k 



Gx — Gy X Gz 



(15) 



Bo ' " |Bo X k| 
Also, we introduce the following wavenumbers 
fcp = A*(7 " l)po^p{RcsToy^ , kr = p.{-y - l)flT{Rcs)~^ , 
'^ify = [m(7 - '^)Kl\^'^{RCsPo), 

kii^ = [p(7 - l)Ki]-\RcsPo). (16) 

Now, we can write the dispersion equation using the follow- 
ing non-dimensional quantities, 

kr k k 

V 



r = — (Jp = 



Also, possible effects of CRs diffusion are appeared in 
our final dispersion relation through the non-dimensional 
parameters (j) a^nd tp as 

, _ Per-O 

po 



(17) 



/ / ^ 

rvn 



(18) 



where Vo = (icr ~ i)kp{cs/ k^^)^^ . Defining a new wavenum- 
ber fcc — (7cr — 1)~^{cs/k\]), then we have tpo = kp/kc- 
Also, this parameter can be re-written as i/)o ~ f^i'Jcr — 
l)ft:||/CsTcooi, where TcooI is the cooling time-scale. Obviously, 
the non-dimensional parameter i/jq shows the level of diffu- 
sion of CRs. When there is no diffusion and the thermal 
gas and CRs are well coupled together, we have i/) = 0. But 
as this parameter increases CRs are diffusing more along 
the magnetic field lines, the level of the coupling becomes 
weaker. In our analysis, these non-dimensional parameters 
are the key input parameters to explore possible effects of 
CRs on the thermal instability. 

Therefore, the characteristic equation becomes 

+ (o-T +(TK + IpO^'" + blKicTT + CTif ) + 1 + a 
+lcr<j}/'y]T'^ + [(1 -I- 7a + 7cr<^)(o-T + O k) h ~ ^pH 

-f i/'C(i + ")]r^ + [V'C(i + iq)((jt + okMi - V'Cc^p/t 

+ qC(1 + 7cr<^/7)]r^ + K(l + '~iar^)(pT + O k) 1 1 

~ aCGph + a^CW + ctipCi'^T + <TK ~ <Tp)h = 0, (19) 

and C, — cos^ 9 and 9 is the angle between Bo and k. Also, 
we have cfk = o-Ky C + (^k^ (1 — C) B.nA a = (wa/cs)^, where 
VA is the Alfven velocity. 



4 ANALYSIS 

Equation ((19} describes magneto-thermal waves and un- 
stable modes in a thermal system with CRs. If we set 
= t/j = 0, this algebraic equation reduces to the stan- 
dard characteristic equation of magneto-thermal instability 
without CRs (e.g.. Field 1965). We are mainly interested in 
roots of equation (|19p that are real and positive (i.e., con- 
densation modes). Qualitative analysis of the roots of the 
algebraic polynomial equations is generally done by Hur- 
witz analysis. In a case without magnetic field. lWagner et al] 
(2005) derived conditions of the stability of the system when 
3> 1 or <C 1. However, we found it difficult to do such 
an analysis analytically and present the stability criteria in 
a simplified closed form (except for = 0) because of the 
complicated coefficients of the characteristic equation (|19p . 
B ut we study domains of instability numerically. In analogy 
to lFieldl(|l965l ). we show domains of the instability based on 
the characteristic equation (|19p in a system of coordinates 
defined by Op and gt- Also, in order to study the influence of 
various parameters on the growth rate of the condensation 
modes, we solve the dispersion equation (|19p numerically to 
locate the roots of the equation against the wave number for 
several values of the parameters. 

Domains of the stability are shown in Fig. [T] For sim- 
plicity we consider cases without thermal conduction. In Fig. 
[1] the vertical axis is related to the density derivative of the 
net cooling function, i.e. (jp, and the horizontal axis is the 
normalized form of the temperature derivative of the cool- 
ing function, i.e. ot- The nondimensioanl parameter of the 
diffusion of CRs is i/) = 100 and the ratio 4> varies as 0.1, 0.4 
and 0.9 and also we assume a = 1, 7 = 5/3 and 7c,- = 4/3. 
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Figure 2. Growth rate of thermal condensation mode versus 
wavenumber of the perturbations when a = 1, 7 = 5/3, 7cr = 
4/3, tpo = 0.1, <Jt/<^p = 1/2 and crpan = 0.01. Dashed line repre- 
sents growth rate in the case without CRs. Each curve is labeled 
by 0. 



Figure 3. Growth rate of thermal condensation mode versus 
wavenumber of the perturbations when a = 1, 7 = 5/3, fcr = 
4/3, </> = 0.6, (7t/<^p = 1/2 and rjpaK = 0.01. Each curve is 
labeled by ipo- 
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Top, middle and bottom plots of Fig.[l]are corresponding to 
^ = 0.1, 0.01 and 1, respectively. Light gray regions in this 
Figure are domains of the stability. In order to make easier 
comparison, the area between two dashed lines corresponds 
to the domain of the stability in a case without CRs. We 
see that domains of the stability extend due to the existence 
of CRs, in particular, for nearly transverse perturbations 
this extension is very significant. There are stable regions 
for nearly perpendicular perturbations, for which parallel 
perturbations are not stable. But in the absence of CRs, 
domains of stability are the same for all the cases. There- 
fore, existence of CRs may lead to the formation of highly 
elongated or even filamentary structures. 

As we mentioned for transverse perturbations, i.e. = 
0, we can present a simple analytical criteria for stability 
of the system. In this case, the characteristic equation p9p 
reduces to a third degree algebraic equation for which one 
can simply do Hurwitz analysis. Criteria of the stability is 
the following inequality 



1 + 7cr0 + 



(20) 



Obviously, in the absence of CRs and magnetic field the 
above c ondition red uces to the classical thermal stability 
criteria (|Fieldlll965h . 

Now, we study condensation modes numerically. We 
take parameters a = 1, ar l<Jp = 1/2 and GpOK = 0.01 for 
easier comparison to iFieldl (| 19651 ). Also, the other parame- 
ters are 7 = 5/3 and 7cr = 4/3. Effects of the initial ratio of 
the CRs pressure to the gas pressure on the growth rate of 
condensation modes is explored in Fig. [2]by fixing t/jq = 0.1 
and changing (j> from 0.1 to 0.6. In the case of no CRs, the 
growth rate is represented by a dashed curve. Top plot is for 
perturbations that are nearly perpendicular to the initial 
magnetic field line, i.e. C, — 0.01. But the bottom plot shows 
growth rate for the cases with the perturbations along the 
initial magnetic field line, i.e. C = 1- Note that each curve is 
labelled by corresponding ratio (j>. Reduced growth rates are 
seen in Fig.[2]due to the existence of CRs, irrespective of the 
direction of the perturbations. Also, reduction of the growth 
rate is larger as the ratio of the unperturbed CRs pressure to 
gas pressure increases. This stabilizing effect of CRs can be 
easily understood by noting to the fact that CRs pressure is 
acting as an extra support to help the gas pressure against 
the unstable perturbations in the momentum equation. 

Comparing to a case without CRs, for a fixed ratio (j>, 
the reduction to the growth rate of the nearly transverse 
perturbations is more than a case with perturbations along 
the magnetic field line according to Fig. [2] In our analysis, 
CRs particles are diffusing along the magnetic field lines. 
Thus, diffusion of CRs is more or less negligible for nearly 
transverse perturbations which implies a maximum pressure 
support by CRs in the momentum equation. 

The profile of the growth rate of condensation mode 
reaches to its largest value for a wavenumber fcmax which 
depends on the input parameters. Fig. [2] shows that in the 
absence of CRs as the direction of the perturbations tends 
to be along the magnetic field line, the maximum wavenum- 
ber fcniax decreases. But when dynamics of CRs is included, 
the wavenumber fcmax depends on the ratio <f> and the non- 
dimensional diffusion coefficient ■0o. For nearly transverse 



perturbations. Fig. [2] shows that the wavenumber fcmax de- 
creases when the ratio (p increases. But this behavior is re- 
versed when perturbations are along the magnetic field lines, 
i.e. fcniax increases with the ratio (j>- 

Fig. [3] shows the effect of CRs diffusion along the mag- 
netic field lines. All the input parameter are the same as 
Fig. [2I except for ipo which varies from to 100 and the 
ratio <j> is kept to be fixed, i.e. (j) ~ 0.6. Growth rate of 
the condensation mode increases with the diffusion of CRs. 
This destabilizing effect of diffusion is enhanced for the per- 
turbations along the magnetic field lines. For this type of 
perturbation the maximum wavenumber fcmax is decreasing 
as more CRs particles are diffusing. Note that when diffu- 
sion of CRs is efficient, the growth rate tends to the profile 
of a case without CRs. This behavior is more evident for the 
perturbations along the magnetic field lines when ^/jq is large 
according to Fig. [3] In other words, compression along the 
magnetic field lines due to the magnetothermal instability 
is stronger when diffusion of CRs are considered. 



5 DISCUSSION 

Interstellar CRs can impinge on the structure formation in 
interstellar medium either by contributing to ionization or 
constituting as an extra source of heating or even through 
dynamical coupling to the plasma. In this study, we ana- 
lyzed dynamical effects of CRs on the unstable modes in 
magnetothermal systems. CRs and the plasma are consid- 
ered as two different fluids and our linear analysis implies a 
stabilizing effect due to the existence of CRs. Since CRs are 
diffusing along the magnetic field line, the level of stabiliza- 
tion of the condensation modes decreases, in particular for 
the perturbations along the magnetic field lines. In fact, in 
the classical analysis of magnetothermal instability one can 
show that a purely tr ansverse m agnetic field can prevent 
thermal condensation (jFieldl l 1965*1 . On the other hand, CRs 
can diffuse along the magnetic field line and not perpendicu- 
lar to it. Therefore, the magnetic and CRs pressures add up 
to the thermal pressure for transverse perturbations where 
diffusion of CRs is negligible. This implies a more magne- 
tothermally stable system. However, CRs pressure reduces 
for the perturbations along the magnetic field line because 
of the diffusion of CRs particles. Since the stabilization ef- 
fect of CRs is stronger for the transverse perturbations, we 
may expect formation of elongated clouds via thermal in- 
stability in the presence of CRs. Effectiveness of CRs on the 
reduction of the growth rate is anisotropic and CRs diffusion 
can be a key factor in the final alignment of the elongated 
clouds. 

We can also compare our analysis with similar previous 
studies like KP83. Although the main goal of KP83 is about 
possible role of CRs in ISM, in comparison to our study there 
are significant differences that we summarize them here: 
(a) The present paper studies thermal instability, and so, 
the net cooling function is considered. But cooling function 
is completely neglected in KP83 and the same for the ther- 
mal conduction. Adiabatic variations of the gas component 
is considered in KP83 instead of the complete form of the 
energy equation (see, eqs. 22 and 23 of KP83). That is a 
significant difference. Because in our analysis thermal insta- 
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bility is occurring on cooling time scale, but in KP83 there 
is not thermal effects due to the cooling of the system. So, 
KP83 can not address thermal condensation modes with CRs 
in a typical ISM (such as HI regions) which may have a vital 
role in theory of star formation. But our analysis discusses 
about thermal condensation modes including dynamical role 
of CRs. Therefore, KP83 and the current paper are address- 
ing somewhat different systems, but leading to a consistent 
picture for the role of CRs. KP83 showed that CRs oper- 
ate as an stabilizing factor in Parker instability. Our paper 
shows that CRs play the same role but within the context 
of thermal instability in ISM even at scales comparable to 
HI regions. 

(b) In the light of the above point, our initial state is also 
different. While we start with a homogenous initial config- 
uration, KP83 starts with an initial configuration in which 
the gas is supported in the vertical direction by both grav- 
ity and magnetic and thermal pressures. Then, their initial 
configuration depends on the vertical spatial coordinate (see 
Figure 1 of KP83). This initial set up is similar to the classi- 
cal Parker analysis. In our analysis, self-gravity is neglected. 
Because we are interested in formation of structures due to 
thermal instability in systems where self-gravity has a neg- 
ligible role (like HI regions). 

(c) Our initial state is actually corresponding to a thermal 
equilibrium configuration for which the net cooling function 
is zero. But initial state of KP83 is defined by the balance 
of different forces in the vertical direction. 

(d) Considering the above points, we determined domains of 
the stability in Fig. [T] numerically (see also equation (|20p '). 
Since we have considered the complete form of the energy 
equation including the net cooling function, the stability re- 
gions are mainly determined by the net cooling function and 
the parameters relating to the CRs. We show that when the 
net cooling of the system is included the stability criteria is 
not a trivial issue. But since KP83 did not take into account 
the energy equation for the gas component, their stability 
criteria has been written in terms of a critical polytropic ex- 
ponent. Therefore, one can not use their stability condition 
for a system with a given net cooling function as we have in 
our study. 

We showed that the level of diffusion of CRs particles is 
determined by the non-dimensional coefficient t/jq- But de- 
pendence of this parameter on the physical parameters of the 
system is a complicated function of the net cooling function. 
Having a fixe d value for the diffusion coefficient k \\ around 
10^* cm^ s"^ l|Berezinskn et al.lll990l : |Ptuskinll200 j) . our key 
parameter ^/jq depends on the physical variables that control 
the rate of cooling such as temperature and ionization frac- 
tion. For simplified cooling systems, we can derive depen- 
dence of -00 on the param e ters o f the system. For instance, 
ISchwarz. McCrav fc SteinI (|l972l ) studied formation forma- 
tion of clouds via thermal instability in ISM. They showed 
that the cooling time scale of a low-density plasma of cosmic 
abundance can be approximated as TcooI ~ T/nx for tem- 
peratures within the range of 100 K to 10* K. Here, n is the 
number density of particles and x is the fraction density of 
the electrons, i.e. x — Ue/n. Then, we have ijjo oc Kypx/T^. 
Therefore, in cold systems, one may expect an efficient dif- 
fusion of CRs along the magnetic field lines which implies 



negligible dynamical effect of CRs on the condensation mode 
(see Fig. [3]). But as the temperature of the system increases, 
the level of diffusion decreases. Note that K|| is assumed to 
be fixed and independent of the properties of the system in 
the above argument. However, there are some points that 
make the problem more complicated. An unstable magne- 
tothermal system may become turbulent at nonlinear regime 
of the evolution. Then, diffusion of CRs particle may be af- 
fected by the level of turbulence within the system. 

In the light of our results, we think, the formation of 
structures in the protogalactic halo environment can be re- 
analyzed by considering dynamical role of CRs. For exam- 
ple, B^k et aL (2006;) found that dense clumps first form 
out of hot background gas by thermal instability in the pro- 
togalactic environment. For such systems, the cooling time 
scale and the sound spee d are estimated a s t^ooI = 2 x lO'^ 
yr and Cs = 198 km s"^ (|Baek et all 120061 ). Then, the non- 
dimensional parameter diffusion parameter becomes Vo ~ 
0.04. This implies that the diffusion of CRs along the mag- 
netic field line is going to be negligible and the growth rate 
of condensation mode is significantly reduced depending on 
the ratio parameter 0. Therefore, the rate of the structure 
formation in the protogalactic halo environment via thermal 
instability is reduced because of the existence of CRs. This 
encouraging result can be studied further via numerical sim- 
ulations of thermal instability with CRs for the protogalactic 
systems. 

In another study, iKim. Kim fc Ostrikeij l|2008t ) studied 
galactic spiral shock waves with thermal instability. They 
showed that initially uniform gas rapidly separates into 
warm and cold phases as a result of thermal instability and 
also for ms a quasi - stead y shock that prompts phase tran- 
sitions (|Kim et al.l |2008| ). Although existence of CRs and 
th eir stabilizing effe ct may not change qualitative scenario 
of iKim et al.] (|2008| ). we can expect the initial gaseous disc 
evolves into warm and cold phases with much slower rate 
comparing to a similar system without CRs. This is an in- 
teresting research topic for future. 



ACKNOWLEDGMENTS 

I thank the anonymous referee for comments leading to a 
greater clarity of the paper. I am grateful to Luke Drury and 
Tom Hartquist for their useful comments and suggestions. 



REFERENCES 

Back C. H., Ryu D., Kang H., Kim J., 2006, ApJ, 643, L83 
Begelman M. C, Zweibel E. C, 1994, ApJ, 431, 689 
Berezinskii V. S., Bulanov S. V., Dogiel V. A., Ptuskin 
V. S., 1990, Astrophysics of cosmic rays. Amsterdam: 
North-Holland, 1990, edited by Ginzburg, V.L. 
Drury L. O., 1983, Reports on Progress in Physics, 46, 973 
Drury L. O., Voelk J. H., 1981, ApJ, 248, 344 
Elmegreen B. G, Elmegreen D. M., 1986, ApJ, 311, 554 
Fatuzzo M., Adams F. C, Melia F., 2006, ApJ, 653, L49 
Ferriere K. M., 2001, Reviews of Modern Physics, 73, 1031 
Field G. B., 1965, ApJ, 142, 531 



© 0000 RAS, MNRAS 000, 000-000 



8 M. Shadmehri 



Field G. B., Goldsmith D. W., Habing H. J., 1969, ApJ, 
155, L149 

Fukue T., Kamaya H., 2007, ApJ, 669, 363 

Gaisser T. K., 1990, Cosmic rays and particle physics. Cam- 
bridge and New York, Cambridge University Press, 1990, 
292 p. 

Goldsmith D. W., Habing H. J., Field G. B., 1969, ApJ, 
158, 173 

Hanasz M., 1997, A&A, 327, 813 
Hanasz M., Lesch H., 2000, ApJ, 543, 235 
Kim C.-G., Kim W.-T., Ostriker E. C, 2008, ApJ, 681, 
1148 

Ko C.-M., Lo Y.-Y., 2009, ApJ, 691, 1587 

Kuwabara T., Ko C.-M., 2006, ApJ, 636, 290 
Kuwabara T., Nakamura K., Ko C. M., 2004, ApJ, 607, 
828 

Kuznetsov V. D., Ptuskin V. S., 1983, A&SS, 94, 5 

Lo Y.-Y., Ko C.-M., 2007, A&A, 469, 829 

Mac Low M.-M., Klessen R. S., 2004, Reviews of Modern 

Physics, 76, 125 
Mouschovias T. C, 1974, ApJ, 192, 37 
Parker E. N., 1966, ApJ, 145, 811 
Ptuskin V. S., 2001, Space Science Reviews, 99, 281 
Scalo J., Elmcgrccn B. G., 2004, ARA&A, 42, 275 
Schlickeiser R., 1989, ApJ, 336, 243 
Schlickeiser R., Lerche I., 1985, A&A, 151, 151 
Schwarz J., McCray R., Stein R. F., 1972, ApJ, 175, 673 
Wagner A. Y., Falle S. A. E. G., Hartquist T. W., Pittard 

J. M., 2005, A&A, 430, 567 
Yusef-Zadeh F., Wardle M., Roy S., 2007, ApJ, 665, L123 



© 0000 RAS, MNRAS 000, 000-000 



